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Abstract 

We calculate next-to-leading QCD corrections to the decay H + — > ud for 
generic up and down quarks in the final state. A recently developed algorithm 
for evaluation of massive two-loop Feynman diagrams is employed to calculate 
renormalization constants of the charged Higgs boson. The origin and summa- 
tion of large logarithmic corrections to the decay rate of the top quark into a 
lighter charged Higgs boson is also explained. 
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1 Introduction 



The accelerators planned to be built in the near future will provide an insight into 
physics at TeV energy scale and thus probe a region especially interesting from the 
point of view of electroweak interactions. Therefore we observe recently an increased 
interest in various aspects of phenomenology of the electroweak symmetry breaking, 
in the framework of the Standard Model and its extensions, which predict one or 
more doublets of Higgs bosons. 

As far as the experimental detection of the Higgs particles is concerned, detailed 
knowledge of their decay properties is of special interest. In the present paper we 
concentrate on hadronic decays of charged Higgs bosons, predicted e.g. by the Min- 
imal Supersymmetric Model (see ref. |3] for a review and further references). In 
the case of the Standard Model Higgs boson hadronic decays have been analyzed in 
a number of publications. QCD and QED corrections were first calculated in 0, 
where it was noted that as the ratio of mass of the decaying scalar particle to that of 
fermions in the final state increases, the one-loop corrections diverge logarithmically. 
This problem was solved by renormalizing mass of quarks at the energy scale equal 
to the mass of the decaying particle and thus absorbing the large correction into the 
tree level decay rate expressed in terms of the running mass. This approach was 
subsequently extended by means of the renormalization group technique in ref. |3|. ffl 
(see also ref. |J for a clear summary), where next-to-leading QCD corrections were 
summed up in the case of large mass of the Higgs boson. The effect of three-loop 
QCD corrections was first calculated in ref. |J and further analyzed in ref. ||. Lead- 
ing logarithmic approximation and various ways to parametrize the next-to-leading 
corrections to the Standard Model Higgs boson are subject of several recent studies 

In the minimal extension of the Standard Model suggested by supersymmetry 
there are two Higgs boson doublets, and one of the charged modes becomes physical. 
It is therefore of great importance to look at the phenomenology of such charged scalar 
particles, since their discovery would give information about the theory underlying the 
Standard Model. Various radiative corrections to hadronic decays of charged scalar 
particles have been calculated and published recently. One-loop QCD corrections and 
the sum of leading logarithms can be found in [ffDL [TTJ] and the leading electroweak 



effects in the limit of large mass of the top quark in [0, [U| (see also [|14|j). 

The purpose of the present paper is to apply the renormalization group technique 
to calculate next-to-leading QCD corrections to the process H + — > ud, where u and d 
represent generic up and down type quarks respectively. In section g we describe the 
theoretical framework of this calculation, and in [3] we present the calculation of the 
two-loop mass and wave function renormalization constants of a charged Higgs boson. 
The results for the corrected rate of the Higgs decay are presented in section f|. In 
section ^| we digress to discuss an application of the summing of leading logarithms 
to the closely related process t — > H + b. Our conclusions are given in section 0. 
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2 Decay rate and operator product expansion 



As has been demonstrated in refs. || £§], the hadronic decay rate of a Higgs boson 
can be represented by Q 



r(# H 



ud) = T H — Im C (q 



where Co is the coefficient of the unit operator in the operator product expansion of 
the correlator function of the scalar currents [], 

i / d D xe 1 ^ T[J H (x)J H (0)] = £C<(g 2 )0<, (2) 
J d,i 

and the scalar current in the present case is defined by 

J H = Z[ — u {aR + bL) d. (3) 

In this formula Z[ denotes the renormalization constant of the charged Higgs-fermion 
vertex, R and L are the chiral projection operators, R = (l+j 5 )/2 and L = (1— 7s)/2, 
and the coefficients a and b depend on the specific model. We will consider two models 
characterized by the absence of flavour changing neutral currents, described in detail 
in where references to original papers can also be found. In model I we have 

a = v2 m u cot j3, b = — Vz m d cot (3, (4) 

whereas in model II, which corresponds to the Higgs sector of the Minimal Super- 
symmetric Standard Model: 

a = V% m u cot (3, b = y/2 m d tan f3, (5) 

and tan (3 is the ratio of vacuum expectation values of the two Higgs doublets. 

We now describe the procedure of calculating next-to-leading QCD contributions 
to the coefficient function Co(q 2 ). Following ref. || and using methods described 
in [H| one derives the renormalization group equation for the space-like values of 
the argument of Co {q 2 < 0). The solution to this equation is found in form of an 
expansion in the QCD coupling constant g s and in the ratios of masses m u ^/mH 
(below we shall also use a common notation for the quark masses, m q = m U) d). 
Keeping only first two terms in the mass expansion is justified in the region far above 
the threshold of the ud production; near the threshold not only this expansion is 
insufficient, but also the perturbative treatment of the vacuum polarization diagrams 
is questionable. 

While the general method follows closely ref. and needs no further discussion 
here, we want to concentrate on the novel feature of the present calculation, namely 
on the computation of the renormalization constants of the charged Higgs boson. 

3 We shall omit "plus" in notations related to H + ; for example, itih = rriH+, Th = r#+, etc. 

4 The sum on the r.h.s. of (2) goes over the indices d and I, where d denotes the canonical 
dimension of a local operator O l d while I labels independent operators with the same d; there is only 
one (unit) operator with d = 0: Oq = 1. 
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3 Derivation of renormalization constants 



We work in D = 4 — u dimensional space, considering 75 to be anticommuting with 
other 7- matrices, 75 = 1. We may use such a scheme because there is no anomaly 
problem in the case considered (see, e.g., ref. |[15|| ). Solution of the renormalization 
group equation requires knowledge of 1/u poles of the quantities 

S = g- 2 (1-Z' 3 ) and T — g~ 2 (l — Z mH Z' 3 ), (6) 

with Z' 3 and Z mH denoting wave function and mass renormalization constants of the 
charged Higgs field. Hence we have to calculate divergent parts of the self energy 
diagrams as depicted in figure [|. It has to be noted that in the present case of two 
different masses of quarks in the loop, diagrams corresponding to figs. |l](b) and |I](d) 
should be considered together with their counterparts with corrections on the other 
quark line. Their sum is of course symmetric under m u <-> m^. 

A modification of the method developed in ref. |T7| enabled us to obtain exact 



expressions for divergent parts of all diagrams presented in figure [I], that are valid for 
any values of the external momentum k, m u and m d . Since in the present paper we 
are mainly interested in the expansion in m 2 q /m 2 H , we expand the relevant integrals 
in m 2 /k 2 , keeping the k 2 and the k 2 (m 2 /k 2 ) = m 2 terms only. In the one-loop order 
(see figure |l|(a)) we get 



ig 2 N c 



c 



(ATiY2m 2 w u 



{Aabm u m d + (a 2 + b 2 ) (2(m 2 + m 2 d ) - k 2 ) } , (7) 



while the sum of all two-loop-order contributions (see figures |l|(b,c,d,e) together with 
counterparts) yields 



(An)Hm 2 w uj 



[a 2 + b 2 



™ (4(m 2 + m 2 ) - k 2 ) + (5A; 2 - 8(m 2 + m 2 ; 
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+16abm u m ( i ( 1 



In the above formulae the colour factors are Afc = 3 and Cp — 4/3. It is remark- 
able that the terms containing ln(m 2 ) and ln(— k 2 ) (occurring in the expressions for 
separate diagrams) disappear in the whole sum @. In particular, this fact enables 
us to consider analytic continuation to time-like values of the momentum without 
difficulty It should be also noted that the factors of 7r _a; / 2 r(i + |o;) are included into 
the definition of coupling constants g 2 and g 2 , as it is usually done in the framework 
of the MS scheme (this is also equivalent to a re-definition of 1/u poles). 

We can now calculate coefficients S n and T n of l/u n poles in S and T including 
terms of the order of g 2 : 

T _ N c 2abm u m d + (a 2 + b 2 )(m 2 u + mj) ( g 2 M) \ 
1 ~ (4vr)2 I (4tt) 2 * ' 
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N c a 2 + b 2 4g 2 

2 " (An) 2 m 2 w (4tt) 2 ' { ' 

N c 2abm u m d + (a 2 + b 2 )(m 2 u + m 2 d ) \Qg 2 s 
2 (4tt) 2 m 2 w m 2 H (4tt) 2 ' 1 ' 

with s^ 1 ' = y an d t = f • We have displayed S 2 and T 2 because they illustrate nice 
agreement of our calculation with equations following from the renormalization group 
analysis 0, |TB[. Namely, both these quantities can be found in the lowest relevant 
order of perturbation theory from Q 

Si = j£y 2 7<? Si + 0{g% T 2 = 2-jL 7 g) T a + 0(^ 4 ), (13) 

and we reproduce our expressions by putting 7$^ = —8. 



4 Corrected decay width 

Here we shall present the formula for the decay rate of the charged Higgs boson 
including next-to- leading corrections. We define C = In (^#/Aq CD ) and obtain: 



H 



N c g 2 m H 

327Tmfy 



£7^ /A) 



X 



a 2 + P 



1 + A7glln£ + 1 



1% 



C 



(3 C 



where the mass correction 5 is 



Im 



+ 5 



(14) 



— i ~ar \^bm u m d + (a 2 + 6 2 )(m 2 + m 2 



- (2abm u m d + (a 2 + b 2 )(m 2 u + m 2 , 

2 / 



x 



Pi 



c 



p £ 



t (1) + 4- 



1 

00 



Im 



(15) 



We have expressed the decay width of the charged Higgs boson in terms of the renor- 
malization group invariant masses of the quarks m q (q = u,d), and the coefficients a 
and b defined by equations (f|) and @ with m q replaced by m q .f\ Although eqs. ( |14"D 



5 Here and below, 7^ and (3 n correspond to the coefficients of expansion (in g s ) of the anomalous 
dimension of mass, 7m(<7s), and the beta function, j3(g s ). In the normalization used (see, e.g., ref. |3| 



102-fiV F , 7 L 0) 



%7- 



and references therein) we have: (3q = 11 — ^Np, (3\ 
where Nf is the number of quark flavours. 

6 If we formally put rh u = frid = a = b = 
partial decay width of a neutral Higgs into qq pair (see eqs. (3.21)-(3.22) in M). 



(i) 



4U4 I AT 

- — - + —N F , 



in eqs. (|lj) and Jl5|), we get the correct result for 



4 



and ( pT5[) correspond to partial decay width H + — > ud, it is clear that the main con- 
tribution to the sum over generations will be given by the term with maximal quark 
masses allowed. 

Following the method of ref. |]3], |18| we use the threshold condition stating that the 



running mass of the quark at the energy scale of production of a pair quark-antiquark 
is equal to half this energy; from this condition we obtain: 



2 _ ™-2 ^""q 



m„ = m„ In 



4m 2 x 7 '° >///3 ° 



q q \ A 2 / 

f A 7 Wlnln(4^/A^ CD ) 1 /ft A 1 1 

1 + /3 3 ln(4ro*/A^ CD ) + ft 7m J m(4m 2 /A 2 QCD ) J ' 1 } 

We now want to visualize the magnitude of the leading and next-to-leading order 
corrections. The leading order correction can be obtained from equations (|TJ]) and 
( |15D by dropping all the terms divided by C = In (m 2 H / Aq CD ) . In the case of model 
II, our formulae in the leading logarithmic approximation reproduce the results of 



ref. |L0| , for both reactions H + — > cs and H + — > £&. To our knowledge, the corrections 
in model I have not been analyzed so far even in the leading order. We present both 
the leading and the next-to-leading logarithmic corrections to the rate of the decay 
H + — > tb in fig. 0(a-c) for the value of tan/3 = 2, for which this decay has similar rate 
in both models (we take m t = 150 GeV and m& = 4.5 GeV). It turns out that although 
the leading corrections decrease the expected width of the charged Higgs boson, the 
next-to-leading terms can increase it by a large factor, especially for the light mass 
of the Higgs. On the other hand, this effect may be an artifact of the choice of the 
threshold condition. The sensitivity of the result to this choice is shown in fig. |2|, 
where we present the results using two different conditions: (a,b) m(4m 2 ) = m and 
(c) m(m 2 ) = m (where m denotes the running mass of the quark). The dependence 
on the initial condition has also been discussed in much detail in ref. 

Size of corrections also strongly depends on the value of tan /3, as can be seen 
in fig. |3|. For this plot we have used the usual condition m(4m 2 ) = m, and we see 
that for large tan/5, where the decay is dominated by the coupling proportional to 
the bottom quark mass, next-to-leading corrections slightly decrease the effect of the 
leading ones. Therefore, the effect of increasing the width mentioned above is due to 
the corrections to the mass of the top quark; this may be a signal of insufficiency of 
the expansion in m q jmn for lighter Higgs bosons. We are going to address this issue 
in future (see also ref. 



5 On the decay t — > H + b 

The same term in the Lagrangian which is responsible for the decay of the charged 
Higgs boson into quarks will also enable a sufficiently heavy top quark to decay into a 
bottom quark and an H + . Both electroweak EO, ETJ and QCD EI3, E3, E3] corrections 



to this decay have been studied at the one-loop level. It has been found [24J that in 



the two Higgs doublet model predicted by supersymmetry the relative size of QCD 
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corrections becomes large for growing values of tan (5. We would like to discuss this 
effect here in order to demonstrate that this large correction can be absorbed in the 
Born rate if one uses the running mass of the b quarkQ just like in the case of the 
Standard Model Higgs boson 0, f25f . 

For this purpose we compute the tree level decay rate in the limit of very large 
tan (3 and mass of the top quark: 



G F ml 



v tb \ 2 (i~ x 2 ) 



4 + (l- X 2 )tan/? 



(17) 



where we have introduced the following notation for the ratios of relevant masses: 



— , X 
m t 



m t 



(m H = m H +). 



(18) 



The first order QCD corrections, calculated in ref. p4J, can be expressed in the limit 
of large tan f3 by: 



T m (t^H + b) 



a s G F ml \ V t 



6tt v/2 



7T 



(2G+ - G )tan 2 /5 + 4G_ 



(19) 



and the explicit formulae for the coefficient functions Gi can be found in ref. p4 
Here we only need terms of the order of lne: 



G + 
GL 
Go 



4 



3(1 -X 2 ) lne, 
3 



;i-X 2 ) 2 lne. 



(20) 



Using these expressions we can calculate the asymptotic value of first order corrections 
for large values of tan/3 and m t : 



7T \m 2 t 



(21) 



We see that for a, ?s 0.1, 



4.5 GeV and m t = 100 GeV this correction is of the 



order of —40%, in agreement with diagrams presented in [24]. The size of corrections 
becomes even larger as the mass of the top quark increases, and eventually the one- 
loop corrected rate of decay becomes negative; such large corrections are a sign of 
a breakdown of the perturbation theory. However, it is possible to avoid the large 
corrections by renormalizing the mass of the b quark not on the mass-shell but at the 
energy scale characteristic to the process, which is the mass of the top quark. The 
running mass of the bottom quark at this energy is: 



m b 



In (4m^/A^ CD ) 
In {ml /A 2 



12/(33-2ATp) 



(22) 



qcd; 



7 One of the authors (A.C.) is grateful to the referee of Physical Review D for suggesting this and 
to Sacha Davidson for helpful discussions on this topic. 
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where Np = 6 is the number of quark flavours, and we take Aqcd = 150 MeV (in 
the MS scheme). We can expand the above expression in a series in the coupling 
constant a s , and we find that: 




m b (m 2 t ) « m 6 1 + — In — § . (23) 



It can now be seen from the formula (|21|) that for large tan (3 and rrit the one-loop 
corrected decay rate approaches the Born rate expressed in terms of the running b 
quark mass. In figure f| we show the dependence of decay rate of the top quark on 
tan (5. We note that for large values of tan f3 the QCD corrected rate is not much 
different from the Born rate expressed in terms of the running b quark mass (j22|), 
and that we no longer face the problem of unreasonably large corrections. 

On the other hand, for the small values of tan/3 the Born rate remains approx- 
imately unchanged when expressed in terms of the running b mass. The reason for 
this is that the dominant coupling of the quarks to the charged Higgs is m t cot (3 in 
this region and mass of the b quark does not play any important role. The same 
can be said about the analysis of the top quark decay in the non-supersymmetric 
two Higgs doublet model, and it explains why no large logarithmic corrections were 
found there for any values of tan f3 . 



6 Summary 

We have found leading and next-to-leading order corrections to the decay width of 
the charged Higgs boson in the framework of two models. In the case of the leading 
corrections in the model motivated by supersymmetry we confirmed the previously 



published formulae ||10|| ; the remaining results are new. For a heavy Higgs boson or 
for large values of tan (5 we found that the next-to-leading order corrections sizably 
decrease the effect of the leading order corrections, and increase the final result for 
the decay rate. We have also examined the process t — > H + b, explaining the origin 
of large logarithmic corrections found in a previous paper |24fl . 
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Figure captions 



Figure 1: Types of diagrams contributing to the wave function and mass renormal- 
ization of the charged Higgs boson: one-loop diagram (a), quark propagator (b) and 
vertex (c) counterterms, and reducible (d) and irreducible (e) two-loop diagrams. 



Figure 2: Rate of the charged Higgs decay (a) in Model II and (b,c) in Model I for 
two different threshold conditions (see text). Dashed: Born rate, dotted: leading, 
and solid: next-to-leading corrections. 



Figure 3: Rate of the charged Higgs decay as a function of tan (3. Dashed: Born rate, 
dotted: leading, and solid: next-to-leading corrections. 



Figure 4: Rate of the decay t — > H + b for m t = 150 GeV, m H = 80 GeV, m b = 4.5 
GeV and a s = 0.1: Born rate (long dash), rate including first order QCD corrections 
(short dash), and the improved Born rate (solid line). 



9 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9309245vl 




(d) (e) 



Figure 1: Types of diagrams contributing to the wave function and mass renormal- 
ization of the charged Higgs boson: onedoop diagram (a), quark propagator (b) and 
vertex (c) counterterms, and reducible (d) and irreducible (e) two-loop diagrams. 
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Figure 2: Rate of the charged Higgs decay (a) in Model II and (b,c) in Model I for 
two different threshold conditions (see text). Dashed: Born rate, dotted: leading, 
and solid: next-to-leading corrections. 
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tan [5 

Figure 4: Rate of the decay t — >• H + b for m t = 150 GeV, m H = 80 GeV, m b = 4.5 
GeV and a s = 0.1: Born rate (long dash), rate including first order QCD corrections 
(short dash), and the improved Born rate (solid line). 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9309245vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9309245vl 



This figure "figl-4.png" is available in "png" format from: 
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